Anomalous power law of quantum reversibility for classically regular dynamics 
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The Loschmidt Echo M(t) (defined as the squared overlap of wave packets evolving with two 
slightly different Hamiltonians) is a measure of quantum reversibility. We investigate its behavior for 
classically quasi-integrable systems. A dominant regime emerges where M(t) oc t~ a with a = 3d/2 
depending solely on the dimension d of the system. This power law decay is faster than the result 
oc t~ d for the decay of classical phase space densities. 
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The search for quantum signatures of chaos has pro- 
vided much insight into how classical dynamics manifests 
itself in quantum mechanics [1,2]. The basic question is 
how to determine from a system's quantum properties 
whether the classical limit of its dynamics is chaotic or 
regular. One very successful approach has been to look 
at the spectral statistics, in particular the distribution of 
level spacings [3]. An altogether different approach, ad- 
vocated by Schack and Caves [4] , has been to investigate 
the sensitivity of the quantum dynamics to perturbations 
of the Hamiltonian. This approach goes back to the early 
work of Peres [5] and has attracted new interest recently 
in connection with the study of decoherence and quan- 
tum reversibility [6-12]. 

The basic quantity in this approach is the so-called 
Loschmidt Echo, i.e. the fidelity 



M(t) = \(^ \eMiHt)exp(-iH t)\i; )\ 2 



(1) 



with which a narrow wavepacket "0o can be reconstructed 
by inverting the dynamics after a time t with a perturbed 
Hamiltonian H = H + V [5,6]. (We set h = 1.) The 
fidelity quantifies the sensitivity of the time-reversal op- 
eration to the uncertainty in the Hamiltonian, and thus 
provides for a measure of quantum reversibility. 

To date, most investigations of M{t) focused on clas- 
sically chaotic Hamiltonians H and Hq [6-10]. One no- 
table exception it the original paper by Peres [5], who 
noted that the decay of M(t) is slower in a regular sys- 
tem — but did not quantify it further. We will show in 
this article that in a regular system a dominant regime 
emerges where M(t) has a power law decay oc t~ 3d / 2 , 
with an exponent depending solely on the dimension d of 
the system. This power law decay establishes the higher 
degree of quantum reversibility of regular systems com- 
pared to chaotic ones, where M(t) decays exponentially. 



This trend is as expected from classical reversibility (de- 
fined in terms of the decay of the overlap of classical phase 
space distributions [13]). However, we find that quantum 
mechanics plays a crucial role in regular systems by in- 
ducing a parametrically faster power law decay oc £~ 3< V 2 
than the classical one oc t~ d . 

We consider the generic situation of a regular or quasi- 
integrable Hq and a perturbation potential V that has 
no common integral of motion with Hq. (By regular or 
quasi-integrable we mean systems with a phase space 
dominated by invariant tori.) This condition ensures 
that, classically, the perturbation has a component trans- 
verse to the invariant tori almost everywhere in phase 
space. Our investigation will moreover focus on a regime 
of sufficiently strong perturbation (defined below), where 
one expects a fast decay of the perturbation correlator. 
This regime is to be contrasted with the linear response 
regime considered in Ref. [11]. 

We follow the semiclassical approach of Jalabert and 
Pastawski [6]. We start from a Gaussian wavepacket 
M*o) = (™ 2 )- d/4 exp[*p . ( r ' - T ) - \r' Q - r | 2 /2a 2 ] 
and approximate its time evolution by 

exp(-tHt)Mr) = fdT' J2Kf(v,T' ;t)M^o), (2) 
K?(T,T> ;t) = Cy 2 exp[i5f (r,r^) - mt Ms /2]. (3) 

The semiclassical propagator is expressed as a sum over 
classical trajectories (labelled s) connecting r and r' in 
the time t. For each s, the partial propagator contains 
the action integral Sf (T,T ;t) along s, a Maslov index 
)i s (which will drop out), and the determinant C s of the 
monodromy matrix. Since we consider a narrow initial 
wavepacket, we linearize the action in r' — tq and per- 
form the integration over r' Q . After a stationary phase 
approximation, the semiclassical fidelity reads 



2\d 



M(t) = (47TCT 2 ) 



dr^Xf(r,ro;t)*Xf (r,ro;Oexp(- ( 7 2 |p s -p | 2 ) 



(4) 



with initial momentum p s = —dS s /dro. 

Eqs. (2-4) are equally valid for regular and chaotic Hamiltonians, as long as semiclassics applies. Squaring the 
amplitude in Eq. (4) leads to a double sum over classical paths s and s' and a double integration over coordinates r 
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and r'. Accordingly, M(t) = M^ d \t) + M( nd '(t) splits into diagonal (s = s') and nondiagonal (s ^= s') contributions. 
The diagonal contribution sensitively depends on whether Hq is regular or chaotic. Ref. [6] found that M^ d \t) oc 
exp(— Xt) for chaotic dynamics, with A the Lyapunov exponent. We will show that the decay turns into a power law 
M^(t) cx t~ 3d / 2 for regular dynamics. The nondiagonal contribution, on the contrary, is insensitive to the nature 
of the classical dynamics (set by Hq), provided the perturbation Hamiltonian V has no common integral of motion 
with Ho- Ref. [7] found that M^ Dd ^(t) oc exp(— Ft) for chaotic dynamics, with T given by the golden rule spreading 
width of an eigenstate of Hq over the eigenbasis of H. (This golden rule decay requires that T is larger than the level 
spacing A, but smaller than the bandwidth.) We will see that the same exponential decay of M^ nd \t) holds when H 
is regular, so that M^(t) always dominates in the long time limit. Consequently, the fidelity decays exponentially, 
oc exp[— min(r, X)t] for chaotic systems, while for regular systems the decay is algebraic, oc i~ 3d / 2 , as it is then set by 
the diagonal contribution. The golden rule width T still determines the regime of validity of the power law decay via 
the condition T > A. 

Continuing from Eq. (4), and still following Ref. [6], we write M(t) as 

M(t) = (4ircT 2 ) d J dr J dr' ^ C S C S , exp[iSS s (r, r ; t) - iSS s , (r', r ; t)] cxp(-a 2 |p s - p | 2 - cr 2 \p s , - p | 2 ), (5) 

with SS s (r, r ; t) = S^(r,r ;t) — S^°(r,r ;t). Considering first the diagonal contribution M^ d \t), we set s = s' and 
expand the phase difference as 

6S s (v,r ;t)-6S s (r',r ;t)= [ dt W[q(t)] ■ (q(t ) - q'(t )). (6) 

Jo 

The points q and q' lie on the classical path with q(t) = r, q'(t) = r', and q(0) = q'(0) = r . In a regular system, the 
distance between two initially close points increases linearly with time, |q(t) — q'(i)| — (t/t)\r — r'\. Here we depart 
from the exponential divergence oc exp[A(t — t)] assumed in Ref. [6] for chaotic dynamics. 

The spatial integrations and the sums over classical paths in Eq. (5) lead to the phase averaging 

exp{iSS s - iSS' s ) -» (exp(iSS s - iSS' s )} ~ cxp[-i(((55 s - SS' S ) 2 )]. (7) 

Since V and Hq have no common integral of motion, we may expect a fast decay of the correlations, 

(diV^djVidit')}} - USijS(t — t'). (8) 

One then gets 

M (d \t) = (47R7 2 ) d J dr J dY'Y^Clcx^-lU J\t (tA) 2 |r-rf)exp(-2a 2 |p s -po| 2 ) 

= (Ana 2 ) d J dr + J dr_ ^ C 2 exp(-±Ut v 2 _) cxp(-2 ( 7 2 |p s - p | 2 ). (9) 

The Gaussian integration over r = r — r' ensures that r w r', and hence r + = (r + r')/2 w r. One C s is then 
absorbed by a change of variable from r + to p s , and the Gaussian integral over r_ gives a factor oc t~ d / 2 . Finally, 
setting C s w t~ d as is the case in a regular system, we arrive at 

M (d) (t) ott- 3d/2 , (10) 

which is the central result of this paper. The power law (10) holds once the perturbation is strong enough to induce a 
golden rule spreading of the eigenstates of Hq over the eigenbasis of H (which is the range of validity [6,7] of the above 
semiclassical approach), and under the assumption that the perturbation potential varies rapidly along a classical 
trajectory of H . [We used this assumption to average the complex exponential in Eq. (7).] The decay exponent 3d/2 
is insensitive to the choice (8) of a <5-function force correlator. Even a power-law decaying correlator oc \t — t'\~ a (with 
a > 1) results in the same exponent as in Eq. (10). 

The nondiagonal contribution (s ^ s') to Eq. (5) is the same as in Refs. [6,7]. The phase averaging can be performed 
separately for s and s' and one gets 

(exp[i5S s ]) = exp(-±(5S 2 }) - exp J dt J* dP (V[q(t)]V[q(t')]^j . (11) 
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The point q(t) lies on path s with q(0) = ro and q(t) = r. 
If V and Hq have no common integral of motion, the cor- 
relator of V gives the golden rule decay oc exp(— Tt) re- 
gardless of whether H is chaotic or regular [15]. We con- 
clude that for regular systems, the fidelity is dominated 
by the algebraically decaying diagonal contribution. 

In order to check numerically the analytical result (10), 
we have studied the kicked top Hamiltonian [1] 

H = (n/2r)S y + (K/2S)S 2 Z ]T S(t - nr), (12) 

n 

which describes a vector spin of conserved magnitude S, 
undergoing a free precession around the y-axis, which 
is periodically perturbed (period r) by sudden rotations 
around the z-axis over an angle proportional to S z . Be- 
cause S is conserved, Hq is a one-dimensional Hamil- 
tonian (d = 1), with a two-dimensional classical phase 
space consisting of the sphere of radius S = 1. The 
canonically conjugated variables are ((p,cos9), where 8 
and (f are spherical coordinates. 

The classical limit of the kicked top is given by the 
map [1] 

{Xn+i = z n cos(Kx n ) + y n sm(Kx n ) 
y n +i = -z n sa\{Kx n ) + y n cos(Kx n ) (13) 
z n+l ~2-ni 

in the cartesian coordinates x — sin cos tp, y = 
sin0sin<p, and z — cos 9. Depending on the kicking 
strength K, the classical dynamics is regular, partially 
chaotic, or fully chaotic. We consider a kicking strength 
K = 1.1 for which the dynamics is regular for most of 
phase space. We checked that our results are not sensi- 
tive to the value of K, as long as the dynamics remains 
regular. 

The quantum mechanical time evolution after n peri- 
ods is given by the n-th power of the Floquet operator 

F = cxp[-t(K/2S)S 2 } exp[-i(7r/2)S B ]. (14) 

We perturb the reversed time evolution by a periodic 
rotation of constant angle around the x-axis, slightly de- 
layed with respect to the kicks in H 0} 

V = 4>S x ^ j 8{t-nT-e). (15) 

n 

The corresponding Floquet operator is f = 
exp(— i(j)S x )Fo. We set r = 1 for ease of notation, and 
varied S between 250 and 1000 (both H and H n conserve 
the spin magnitude). We calculated the average decay M 
of M(t = n) = \{ij) \{F't) n F n \ip Q )\ 2 taken over 50 to 200 
initial Gaussian wavepackets Vo OI minimal spreading 
(coherent states). 

In Fig. 1 we show the decay of M for S — 1000 and 
different perturbation strengths (p. For weak perturba- 
tions, the decay of M is exponential, and not Gaussian 
as one would expect from first order perturbation theory 



[5] . The reason why the perturbation operator S x gives 
no first order correction is that for K = 1.1, eigenfunc- 
tions of Fq are almost identical to eigenfunctions of S y , 
so that diagonal matrix elements of V vanish in this ba- 
sis. For weak (f>, the local spectral density of states p(e) 
consists then of a delta function at zero energy plus an 
algebraically decaying tail [14]. Because of the absence of 
a first-order correction, the decay of the fidelity is given 
by the Fourier transform of this tail [10]. We numerically 
obtained a decay p(e) oc (e 2 +7 2 /4)~ 1 with 7 oc (j) 1 - 5 . The 
resulting exponential decay oc exp(— jt) of the fidelity dif- 
fers from the golden rule decay oc exp(— Tt) with r oc 4> 2 . 
As (f> increases, the decay of M turns into the predicted 
power law oc t~ 3//2 , which prevails as soon as one enters 
the golden rule regime, i.e. for T/A « (j) 2 S 3 > 1 [7]. One 
therefore expects the power law decay to appear as S is 
increased at fixed 4>, which is indeed observed in the inset 
to Fig. 1. 

We checked that these results are not sensitive to our 
choice of Hamiltonian, by replacing S x in Eq. (15) with 
S 2 (this is the model used in Rcf. [11]) and also by study- 
ing a kicked rotator as an alternative model to the kicked 
top. These numerical results all give clear confirmation 
of the power law decay (10). 

It is instructive to contrast these results for the decay 
of the overlap of quantum wavefunctions with the decay 
of the overlap of classical phase space distributions, a 
"classical fidelity" problem that has recently been inves- 
tigated [9,11,13]. We assume that the two phase space 
distributions po and p are initially identical and evolve ac- 
cording to the Liouville equation of motion corresponding 
to the classical map (13) for two different Hamiltonians 
Ha and H. We consider regular dynamics and ask for 
the decay of the normalized phase space overlap 

M c (t) = Jdxj dpp (x,p;t) p( X ,p;t)/N p , (16) 

where N p = (j dx / dp p ) x l 2 (j dx / dp p) 1 / 2 . 

We have found above that a factor oc t~ d / 2 in the de- 
cay of the quantum fidelity M(t) oc t~ 3<i / 2 originates from 
the action phase difference and is thus of purely quantum 
origin. One therefore expects a slower classical decay 
M c (t) oc C s oc t~ d . In Fig. 2 we show the decay of the 
averaged M c taken over 10 4 initial points within a nar- 
row volume of phase space a = sin 659 Sip, for K = 1.1 
and <j> = 1.7 • 10~ 4 . The decay is M c oc t -1 , and clearly 
differs from the quantum decay oc i -3 / 2 . 

The power law decay prevails for classically weak per- 
turbations, for which the center of mass of p and po stay 
close together. [This is required by the stationary phase 
condition leading to Eq. (4).] Keeping a fixed, and in- 
creasing the perturbation strength <f>, the invariant tori 
of H start to differ significantly from those of H on the 
resolution scale a, giving a threshold 4> c a. Above 
<j> c , the distance between the center of mass of po an d p 
increases with time oc t and one expects a much faster 
decay M c (t) oc exp[— const x t 2 ] for classical Gaussian 
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phase space distributions [13]. Quantum mechanically, 
a = 1/S (the effective Planck constant) and the thresh- 
old translates into <p c ~ 1/S, coinciding with the upper 
boundary of the golden rule regime. As long as one stays 
in that regime, the perturbation will affect the phase in 
Eq. (7), and result in the anomalous power law decay 
oc t- 3d / 2 . 

In conclusion, our investigations of the Loschmidt Echo 
(1) in the generic regime of classically quasi-integrable 
dynamics show that its decay is dominated by the power 
law M(t) oc t~ a . While from purely classical considera- 
tions one expects an exponent a c = d, we semiclassically 
obtain an anomalous exponent a = 3d/2. This is corrob- 
orated by numerical simulations. The power law decay 
is to be contrasted with the exponential decay found for 
chaotic systems, thereby providing for a novel "quantum 
signature of chaos" . 
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FIG. 1. Decay of M for S = 1000, K = 1.1, and 
10 5 4> = 1.5, 4.5, and 10 (thick solid lines from right to left). 
The crossover from exponential to power-law decay is illus- 
trated by the dotted-dashed line oc exp[— 2.56 ■ 10 -5 t] and 
the dashed line oc t~ A ^ 2 . The dotted line gives the classical 
decay oc t -1 . Inset: Decay of M for K — 1.1, <j> — 10~ 4 , 
and S = 250, 500, and 1000 (solid lines from right to left). 
The dashed and dotted-dashed lines indicate the power law 
oc t~ 3 / 2 and exponential oc exp[— 2- 10~ 4 t] decay, respectively. 
These plots show that the t~ 3 ^ 2 decay is reached either by in- 
creasing the perturbation strength <j> at fixed spin magnitude 
S, or by increasing S at fixed <j>. 
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FIG. 2. Decay of the average overlap (16) of classical phase 
space distributions, for the kicked top with K = 1.1 and 
4> = 1.7 ■ 10~ 4 (solid line). The dotted and dashed lines 
give the classical and quantum power law decays oc t~ x and 
oc t~ 3 / 2 , respectively. 
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